Quantum and thermal effects in dark soliton formation and dynamics in a ID Bose gas 
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We numerically study the imprinting and dynamics of dark solitons in a bosonic atomic gas in a tightly- 
confined one-dimensional harmonic trap both with and without an optical lattice. Quantum and thermal fluctua- 
tions are synthesized within the truncated Wigner approximation in the quasi-condensate description. We track 
the soliton coordinates and calculate position and velocity uncertainties. We find that the phase fluctuations 
lower the classically predicted soliton speed and seed instabilities. Individual runs show interactions of solitons 
with sound waves, splitting and disappearing solitons. 
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Dark solitons, or phase kinks, are commonly studied ex- 
citations in atomic Bose-Einstein condensates (BECs) and in 
nonlinear optics (It], and are reminiscent of robust nonlinear 
solitary waves emerging in numerous other physical systems. 
The main research emphasis has been on classical mean-field 
properties of weakly interacting systems that can be accurate 
modeled by the nonlinear Schrodinger or the Gross-Pitaevskii 
equation (GPE) In BEC experiments dark solitons were 
prepared by imprinting a sharp phase jump on the atomic 
cloud by optical potentials by ultraslow light [4], or by 
merging two BECs (3], and the two-dimensional (2D) and 3D 
soliton dynamics were found to be unstable to perturbations 
perpendicular to the soliton velocity (snake instability). Clas- 
sically solitons can also decay due to sound emission gener- 
ated by an anharmonic trapping potential SSBS- 

In atomic BECs solitons are manifestations of quantum me- 
chanics on a macroscopic (or mesoscopic) scale. Soliton prop- 
erties are sensitive to fluctuations of the BEC wavefunction 
and provide an ideal system to study the emergence of clas- 
sical physics from interacting quantum many-body physics 
ifioll . In recent experiments bosonic atoms were confined in 
tight elongated ID traps that, together with optical lattices, 
significantly enhance quantum fluctuations due to stronger ef- 
fects of interactions (Til [13], rendering classical descriptions 
invalid. In this Letter we study quantum properties of solitons 
by considering the phase imprinting process and the resulting 
nonequilibrium quantum dynamics of phase kinks in such sys- 
tems by including the enhanced quantum and thermal fluctua- 
tions of the atoms within the truncated Wigner approximation 
(TWA)|[i3l[ii[il. In ID the snake instability is suppressed 
and quantum features of solitons are easily distinguishable. 
Individual stochastic realizations of TWA represent possible 
outcomes of single experimental runs revealing jittering os- 
cillatory motion, splitting and disappearing solitons (Fig. [T]), 
while the ensemble averages calculated from TWA produce a 
quantum statistical description of the soliton dynamics over 
its entire trajectory. Numerical tracking of soliton coordinates 
at diflTerent times allows us to evaluate the position and ve- 
locity uncertainties of the phase kinks. We find them rapidly 
growing as a function of quantum and thermal depletion, but, 
surprisingly, observe the average soliton speed being reduced 
due to enhanced fluctuations as a result of the nonlinear de- 



FIG. 1 : (Color online) Atom densities of individual stochastic real- 
izations of soliton dynamics in a combined harmonic trap and optical 
lattice at r = 0. Top: Diff'erent single- shot outcomes for the same 
system obtained from the Wigner function projected to the lowest 
energy band in a lattice. The imprinted phase in the corresponding 
classical case 0^ = 1.3 and A^tot = 240. The solitons disappear or 
split. Bottom: The Wigner density for 0^ = tt, Mot - 900 (left), and 
for two diff'erent stochastic trajectories (both with 0^ = Mot = 440) 
that start to oscillate in opposite directions (center and right). 



pendence of the soliton speed upon its phase distribution. In 
a lattice quantum and thermal fluctuations seed the dynamical 
instabilities of the corresponding classical system-in extreme 
cases the soliton motion is solely generated by quantum fluc- 
tuations. Fluctuations also mitigate the unstable trapping of 
solitons in individual lattice sites by stimulating sound emis- 
sion and decreasing the position uncertainty of the solitons. 

TWA has proved particularly useful for the studies of dissi- 
pative quantum dynamics of bosonic atoms in ID optical lat- 
tices when the atom filling factor in the lattice is not too small. 
For example, it qualitatively produced the experimentally ob- 
served damping rate 1.111 of the center-of-mass motion of the 
atom cloud even for significant ground state depletion fisll . 
Here we show that the quantum statistics of the atoms in TWA 
can be synthesized according to the quasi-condensate descrip- 
tion, analogous to the long-wavelength limit of the Luttinger 
liquid theory, providing an improved representation of TWA- 
phase kinks are particularly sensitive to the resulting enhanced 
phase fluctuations. By varying the eflfective interaction pa- 
rameter we can study a smooth transition of soliton dynamics 
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from a classical to a strongly fluctuating regime. An addi- 
tional advantage of TWA over quantum field-theoretical soli- 
ton descriptions is that it can more easily incorporate excita- 
tions of the system far from the thermal equilibrium as typical 
in soliton imprinting experiments. In nonlinear optics TWA 
was used to analyze solitons in fibers In previous theo- 
retical dark soliton studies in BECs thermal atoms were shown 



to damp the soliton motion 111711 and quantum fluctuations 111811 
fill the soliton core lislll^l or aff'ect the statistics |211. 

We assume the atom dynamics to be restricted to ID by a 
strong radial confinement, so that TWA dynamics follow from 
an ensemble of stochastic fields i//w satisfying [13] 



(1) 



with giD = 2fiajra, where a denotes the 5-- wave scattering 
length and ojr the radial trap frequency. Equation ([T]) formally 
coincides with GPE, but here </^w(x, t) is a stochastic phase- 
space representation of the full field operator. The noise is in- 
cluded in the initial conditions of i//w that synthesize the quan- 
tum statistical correlation functions of the atoms and generate 
the fluctuations of TWA dynamics. As phase kinks represent 
defects sensitive to phase fluctuations, we generate the initial 
state noise for the bosonic field operator ilr(x, 0) carefully us- 
ing a ID quasi-condensate description 1I22I1 by introducing the 
density Sp(x) and phase (f(x) operators 



if/(x, 0) = ypo(x) -h Sp(x) exp(/^(x)), 
Sp(x) = ^Jpo(x) ^ {Spj(x)aj -h Sp*(x)a^j) , 



(2) 



where (pj(x) = Uj(x) -\- Vj(x) and Spj(x) = Uj(x) - Vj(x) are 
given in terms of the solutions to the Bogoliubov equations, 
Uj(x) and Vj(x) Here po = No\i//o(x)\^ and i//o(x) is the 

ground state wavef unction with ^0 particles. 

In TWA the initial conditions are given by Eq. (O, with the 
operators (aj, aj) replaced by complex variables (a*, aj) sam- 
pled from the relevant Wigner distribution-in this case repre- 
senting ideal harmonic oscillators in a thermal bath JTJl, such 
that {a*aj}w = hj + ^. Here ^ results from the Wigner repre- 
sentation that returns symmetrically ordered expectation val- 
ues. We consider the total atom number A^tot to be fixed at 
each stochastic run. The ground state population A^o for a par- 
ticular run may then be defined with reference to the depleted 
population A^nc of the same run 
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FIG. 2: (Color online) Soliton dynamics in a harmonic trap without a 
lattice, (a-d) The Wigner density \4fw{xJ)\^ for individual stochastic 
realizations (that represent possible outcomes of single experimen- 
tal realizations) with the same gioMot and 0^ = 2; (a-c) At T = 
forgot = 50,100,900, corresponding to A^nc/A^tot = 0.2,0.1,0.01, 
respectively. The emitted sound pulses from the imprinting process 
interact more with the soliton at low atom numbers and large quan- 
tum fluctuations; (d) Mot = 900, keT jficx)^ = 22 with A^nc/Mot = 0.3. 
Thermal eff'ects are weaker even with large thermal depletion; (e) 
Soliton velocity at 7 = that would follow from cos((0)/2) for the 
expectation value of the imprinted phase (0) (red/top curve) and the 
expectation value of the velocity obtained from (cos(0/2)> (black 
curve) immediately after the imprinting vs depletion; (f) The en- 
semble averages of the soliton positions {x) (solid lines) and the 
quantum mechanical position uncertainties 6x (shaded regions) for 
Mot = 8000, 440, 50 (curves with decreasing amplitudes) with the 
same nonlinearity gioMot. 0c = 2, and T = 0. At Mot - 8000, 5x 
is negligible. Quantum fluctuations increase 6x and soliton damping, 
and decrease the speed. 



a simple description to sample the ground state mode operator. 
In the rest of the paper, when specifying the depletion we re- 
fer to the expectation value A^nc/Mot- We consider the eff'ects 
of varying A^nc/Mot, either by changing T or g id /Mot whilst 
keeping the nonlinearity constant at gioMot = l^^ficoxh^ 
where ojx denotes the axial trap frequency and Ix = ^|f^Jnuo~x. 
Here ^ id /Mot ^ Tint where the eff'ective ID interaction pa- 
rameter 7int = mgw/fi^n and n is the ID atom density. For 
comparison, in the recent experiment IllH . atoms were con- 
fined to an array of ID tubes with gioMot - '^20fiojxlx and 
Mot - 70 in the central tube. The ratio g id /Mot can be experi- 
mentally managed, e.g., by adjusting the radial trap frequency 
or the scattering length as Mot is varied. 

We numerically study the imprinting process of phase kinks 
in a ID bosonic gas and the resulting TWA quantum dynam- 
ics ([T]). We consider the experimental imprinting method (sll, 
where a soliton is generated by applying an additional con- 
stant 'light- sheet potential', of value to half of the atom 
cloud, for time r, so that the external potential in Eq. ([T]) reads 



Mc = J ^^Yj [('^^''^ " ^^^) (l^iWI^ + IVjWP) + \vj(x)\^] , V = mojlx^/2 -h sEr sm^(7ix/d) -h 0(t - t)e(x)V^ . 



(4) 



(3) 

so that for each run we set A^o = Mot -Mc- Here A^nc fluctuates 
at each realization around its average value Mnc, obtained from 
Eq. O with the substitution {ajf- {a^j^j) + 1/2, and A^o fluc- 
tuates around the average value Mq = Mot -Mc- This provides 



Here s denotes the lattice strength in the units of lattice pho- 
ton recoil energy Er = ff'7i^/2md^ with the lattice spacing 
d. In the corresponding classical case the light sheet im- 
prints a phase jump of cpc = V(pT/fi at x = 0, preparing a 
dark soliton. The imprinted solitons were consequently free 
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to evolve in a harmonic trap (Fig. [2]) and in a combined har- 
monic trap and lattice (Fig. [T]). The corresponding classi- 
cal soliton (governed by GPE) oscillates in a harmonic trap 
(s = 0) at the frequency coxl ^ 1 2] with the initial velocity 
v/c = cos(0c/2), depending on the imprinted phase 0c, where 
c{x) = ^jgix)Nxox\il/{x)\^ Im is the speed of sound. The soliton 
is stationary (dark) for (pc = 7i, with a zero density at the kink. 
Other phase jumps produce moving (grey) solitons, with non- 
vanishing densities rig = n cos^(0c/2) at the phase kink, where 
n is the background density, and |v| c for 0^ ^ 0- Hence, 
the soliton speed and depth can be controlled by r or [ 3] . 

In quantum case, soliton trajectories in TWA fluctuate be- 
tween difl'erent realizations. Individual stochastic realizations 
of |0w^P in a harmonic trap (^ = 0) for difl'erent atom num- 
bers and temperatures are shown in Fig.[2ta-d) that represent 
possible experimental observations (quantum measurements) 
of single runs. The phase imprinting creates sound waves ||3|] 
that more strongly interact with the individual soliton trajec- 
tories at low atom numbers (for a fixed nonlinearity). We also 
show the relative eff'ects of thermal versus quantum depletion. 
We find the eff'ect of even a weak T = quantum depletion 
clearly more significant than thermal depletion (evaluated for 
large atom numbers for which 7 = depleted fraction is neg- 
ligible), in both increased deviations from the classical sinu- 
soidal oscillations and in damping. Fluctuations 'damp' the 
dipolar soliton motion such that the oscillation amplitude is 
increased, but the frequencies fluctuate only slightly around 
the classical value of (jj^I V2. In TWA we can ensemble av- 
erage hundreds of stochastic realizations in order to obtain 
quantum statistical correlations of the soliton dynamics. We 
numerically track the position of the kink at difl'erent times in 
individual realizations and calculate the quantum mechanical 
expectation values for the soliton position (x) and its uncer- 
tainty 6x - ^(x^) - (x)^. The Wigner probabilities can be 
transformed to normally ordered expectation values by sub- 
tracting half-an-atom per phonon mode. In Fig. Of) we show 
(x) and 6x for systems with difl'erent values of ground state 
depletion, clearly demonstrating increasing 6x due to quan- 
tum fluctuations [also shown in Fig. [3tb)]. In Fig. Oa) we 
find that 6x at later times is larger than predicted by the initial 
velocity uncertainty alone. We also fit the individual soliton 
trajectories with the curve x{f) - f(t) exp(-yt), where f(t) is 
an undamped sinusoid and find the quantum expectation value 
for the damping rate y ^ -O.OSlaj^ for 7 = 0, Mot = 900, 
increasing with T up to y ^ -O.OlSoJx at T = 22fioj/kB, cor- 
responding to 30% depleted fraction. At T = 0, y is reduced 
with increasing A^tot, with y ^ -0.023^;;^ at A^tot - 8000. The 
qualitative features of the finite-temperature damping for large 
atom numbers (when quantum fluctuations are negligible at 
T = 0) compares to the system analyzed in a difl'erent regime 
in Ref. Hvl] , where noticeable damping only occurs over many 
oscillation periods. 

We find that the fluctuations of the phase jump across the 
soliton dramatically afl'ect the expectation value and the un- 
certainty of the soliton velocity [Fig.[3e,f)]. Similarly to pre- 




FIG. 3: (Color online) (a) The evolution of the soliton position uncer- 
tainty Sx in a harmonic trap as numerically tracked from the soliton 
coordinates at T = for Mot = 440 (Mc/Mot - 0.02; soHd line) and 
the position uncertainty of an ensemble of sinusoids with the corre- 
sponding initial velocity distribution (dotted line), with Sx quickly 
becoming larger than predicted by the uncertainty in the initial ve- 
locity alone; (b) 6x at four times t < S/oj^ [as marked in (a)] vs de- 
pletion, obtained by varying gio /Mot for constant gioMot- The hori- 
zontal dotted line indicates the soliton core size at the trap center; (c) 
The quantum expectation value of the soliton position (x) (solid line) 
and its uncertainty 6x (dotted line) in a combined harmonic trap and 
optical lattice for 5 = 1, T = 0, Mot = 8000 (Mc/Mot ^ 0.001); (d) 
6x [at times indicated with corresponding symbols in (c)] at T = 0. 
In all cases the soliton positions are calculated before solitons break 
or split in individual trajectories. 



vious studies we find that quantum depletion causes fill- 
ing of the soliton core in the ensemble average of atom den- 
sities. The classical relation between the soliton depth and 
speed suggests the depletion might also increase the soliton 
speed, but we show that, in fact, the opposite is true. The 
quantum expectation value of the phase jump (0) decreases 
due to quantum fluctuations. However, the classical expres- 
sion for the soliton speed | cos(0/2)| has a negative curva- 
ture, so that for a symmetric phase distribution the quan- 
tum expectation value Kcos(0/2))| is always smaller than the 
speed resulting from the expectation value of the phase jump 
I cos((0/2))|, and will decrease as the width of the phase distri- 
bution increases. For instance, if the phase distribution has a 
small width A0, the speed is approximately reduced by the 
factor 1 - A0V32, since the average of cos [(0 ± A0/2)/2] 
approximately yields cos(0/2) [l - A0^/32j. Fig.[2e) shows 
that despite (0) decreasing with increasing ground state de- 
pletion, the average initial speed decreases due to the broader 
phase distribution dominating over the efl'ect of decreasing 
(0). We find that the formula Kcos(0/2))| gives good agree- 
ment with the average initial speed, but underestimates (over- 
estimates) it by a few percent at large (small) atom numbers. 
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Different approximations in the stochastic sampHng of the 
initial quantum statistical correlation functions affect the soli- 
ton dynamics. In the Bogoliubov approximation lll4|] the 
imprinted phase distribution is slightly narrower than in the 
quasi-condensate description. For small depletion the differ- 
ence in dynamics between the two approximations is negligi- 
ble, but for quantum depletion A^nc/Mot - 0-2, it is typically 
about 8%, reaching 30% at later times when Sx becomes large. 

We next consider phase kink dynamics in a combined har- 
monic trap and an optical lattice. In classical soliton dynamics 
a lattice introduces damping and instabilities at low velocities 
JtIIsI]. We consider a lattice with nlx/d = 4, corresponding to 
about 21 occupied sites. Shallow lattices (s = 0.25) show no 
qualitative effects in the soliton dynamics, but this changes at 
about s = 0.5. In Fig. [T] individual stochastic realizations at 
s = 1 display strong effects of quantum fluctuations exhibiting 
rich variation that is not always easily captured by quantum 
mechanical ensemble averaging: fast solitons may disappear, 
or split into multiple solitons at later times, even for small de- 
pletion A^nc/Mot = 0.04. We also calculate Sx (Fig. [3]) that 
becomes very large before the break-up or disappearance. 

Stationary solitons at the harmonic trap minimum in a lat- 
tice are at an unstable equilibrium Jvl], but even at low initial 
velocities first remain trapped in the central site before starting 
to oscillate l'9'i (Fig.IB, as it takes time for the soliton to lose 
energy via sound emission and move into the neighboring site. 
In quantum simulations we find the motion seeded by fluctua- 
tions even for (pc = ^ when the expectation value of the initial 
velocity vanishes, corresponding to large fluctuations for the 
soliton tunneling time out of the central site. For (pc = ^ and 
Mot = 900 at r = the tunneling time and its uncertainty are 
of the order of the trap period. At low initial speeds quantum 
fluctuations may even change the direction of the initial veloc- 
ity, generating large variations between individual trajectories 
and completely dominating the dynamics [Fig. [T] (bottom)]. 
Large soliton uncertainties were also seen in systems of small 
atom numbers in the tight-binding limit lEoll . 

Phase kinks may also become trapped in the outer wells at 
the turning points of their trajectories [Fig.[Tl. The increased 
trapping time in the outer lattice sites accompanies larger fluc- 
tuations in tunneling times, so that Sx increases here for sys- 
tems with smaller depletion. Figure Od) shows that the un- 
certainties in position are initially greatest in the low atom 
number and high temperature cases, due to the broad initial 
velocity distribution, but the situation is reversed after the first 
turning point in the soliton trajectories. Hence the surprising 
result that phase kinks in systems with a small ground state 
depletion exhibit larger position uncertainties here than those 
in systems with a large depletion. The increased sound emis- 
sion is also evidenced by the damping rate. A soliton with 
an imprinted phase jump of cpc = 2.6 has a classical damping 
rate -O.OScOx, which increases to the average quantum value 
y ^ -O.lcox as the number of atoms is reduced to A^tot = 900 
at ^ = 1, r = 0, corresponding to A^ns/Mot = 0.01. Faster 
solitons only experience weak damping (e.g. cpc = 2.0 with 
y = -0.04clJx) with little effect due to similar depletion. 



In conclusion, we showed that solitons can provide an ideal 
dynamical observable to probe underlying quantum fluctua- 
tions. The TWA phase- space description can incorporate a 
very large number of degrees of freedom; the dissipative dy- 
namics of solitons or the entire atom cloud 1IT5I1 emerge from a 
microscopic treatment of the unitary quantum evolution with- 
out any explicit damping terms. The advantage is that the fre- 
quently problematic separation of quantum dynamics to 'sys- 
tem' degrees of freedom and 'environment' [10] is not needed 
and classical physical observables naturally and unambigu- 
ously emerge in the formalism. Here our numerical tracking 
of the soliton coordinates provides outcomes of single- shot 
measurements of soliton trajectories as well as a precise quan- 
tum statistical description of the soliton position and velocity. 

We acknowledge discussions with Y. Castin and financial 
support from EPSRC. 
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